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Abstract 

We argue that space-time properties are not absolute with respect to the used frame of reference as is to be expected 
according to ideas of relativity of space and time properties by Berkley - Leibnitz - Mach - Poincare. From this point of 
view gravitation may manifests itself both as a field in Minkowski space-time and as space-time curvature. If the motion of 
test particles is described by the Thirring Lagrangian, then in the inertial frames of reference, where space-time is pseudo- 
Euclidean, gravitation manifests itself as a field. In reference frames, whose reference body is formed by point masses 
moving under the effect of the field, it appears as Riemannian curvature which in these frames is other than zero. For 
realization of the idea the author bimetric gravitation equations are considered. The spherically - symmetric solution of the 
equations in Minkowski space-time does not lead to the physical singularity in the center. The energy of the gravitational 
field of a point mass is finite. It follows from the properties of the gravitational force that there can exist stable compact 
supermassive configurations of Fermi-gas without an events horizon. 



1 Introduction 

The key reason preventing a correct inclusion of the Ein- 
stein theory of gravitation in the interactions unification is 
that gravity is identified with space-time curvature. It is also 
a cause of such unsolved problems of the theory as an op- 
erational definition of the observable variables, the energy - 
momentum tensor problem and gravity quantization. In the 
present paper starting from [Q, [ |l5| ] we consider a likely 
reason of gravity geometrization. We argue that the grav- 
itation properties are not absolute with respect to the used 
frame of reference. In inertial frames of reference gravita- 
tion can be considered as a field in flat space-time, while 
in so called proper frames of reference it manifests itself as 
space-time curvature. 

The author's gravitation equations which realize this idea 
are considered in details. They do not contradict available 
experimental data. The physical consequences resulting 
from the equations differ very little from the ones in gen- 
eral relativity if the distances from the attracting mass are 
much larger than the Schwarzschild radius r g . However, 
they are completely different at the distances equal to r g 
or less than that. A number of new physical consequences 
follow from the equations. 



The Einstein theory of gravitation demonstrates rela- 
tivity of space-time with respect to distribution of matter. 
However, space-time relativity with respect to measurement 
instruments hitherto has not been realized in physical the- 
ory. An attempt to show that there is also space-time rela- 
tivity to the used reference frames for the first time has been 
undertaken in Q, [Q. 

At present we do not know how the space-time geome- 
try in inertial frames of reference (IFRs) is connected with 
the frames properties. Under the circumstances, we simply 
postulate (according to special relativity) that space - time 
in IFRs is pseudo-Euclidean. Next, we find a space-time 
metric differential form in noninertial frames of reference 
(NIFRs) from the viewpoint of an observer in the NIFR 
who proceeds from the relativity of space and time in the 
Berkley - Leibnitz - Mach - Poincare (BLMP) sense. It is 
shown that there are reasons to believe that side by side 
with generally accepted viewpoint on motion in noniner- 
tial frames of reference as a relative motion there can also 
be another point of view. According to this viewpoint the 
metric differential form ds in the NIFR is completely con- 
ditioned by the properties of the frame being used as is to 
be expected according to the idea of relativity of space and 
time in the BLMP sense. 



2 Primary Principles 

The geometrical properties of space-time can be described 
only by means of measuring instruments. At the same time, 
the description of the properties of measuring instruments, 
strictly speaking, requires knowledge of space-time geom- 
etry. One of the implications of it is that geometrical prop- 
erties of space and time have no experimentally verifiable 
significance by themselves but only within the aggregate 
"geometry + measuring instruments". We got aware of it 
owing to Poincare [|J|. It is a development of the idea by 
Berkley - Leibnitz - Mach about relativity of space-time 
properties which is an alternative to the well known New- 
tonian approach. 

If we proceed from the conception of relativity of space- 
time, we assume that there is no way of quantitative descrip- 
tion of physical phenomena other than attributing them to 
a certain frame of reference which in itself is a physical 
device for space and time measurements. But then the rela- 
tivity of the geometrical properties of space and time men- 
tioned above is nothing else but relativity of space-time ge- 
ometry with respect to the frame of reference being used.Q 

Thus, it should be assumed that the concept of the refer- 
ence frame as a physical object, whose properties are given 
and are independent of the properties of space and time, is 
approximate, and only the aggregate "frame of reference + 
space-time geometry" has a sense. 

1 There is an important difference between a frame of reference (as a 
physical device, '! an d a coordinate system (as way to parameterize points 
of space-time) llq] 



3 The Metric Form ds in NIFR. 

By a noninertial frame of reference we mean the frame, 
whose body of reference is fo/-rm/-ed by the point masses 
moving in the IFR under the effect of a given force field. 

It would be a mistake to identify "a priori" a transition 
from the IFR to the NIFR with the transformation of coor- 
dinates related to the frames. If we act in such a way, we 
already assume that the properties of the space-time in both 
frames are identical. However, for an observer in the NIFR, 
who proceeds from the relativity of space and time in the 
BLMP sense, space-time geometry is not given "a priori" 
and must be ascertained from the analysis of experimental 
data. 

We shall suppose that the reference body (RB) of the 
IFR or NIFR is formed by the identical point masses m p . 
If the observer is at rest in one of the frames, his world line 
will coincide with the world line of some point of the ref- 
erence body. It is obvious to the observer in the IFR that 
the accelerations of the point masses forming the reference 
body are equal to zero. Of course, this fact occurs in rela- 
tivistic sense too. That is, if the differential metric form of 
space-time in the IFR is denoted by dr\ and v a — dx a /drj 
is the 4- velocity vector of the point masses forming the ref- 
erence body, then the absolute derivative of the vector v a is 
equal to zero,i.e. 

Diy a /dj 1 = 0. (1) 

(We mean that an arbitrary coordinate system is used). 

Does this fact occur for an observer in the NIFR ? That 
is, if the differential metric form of space-time in the NIFR 
is denoted by ds, does the 4-velocity vector £ Q = dx a /ds 



of the point masses forming the reference body of this NIFR 
obey the equation 

DQ a Ids = (2) 

? The answer depends on whether space and time are 
absolute in Newtonian sense or they are relative in the BLMP 
sense. 

If space and time are absolute, the point masses of the 
NIFR reference body for an observer in this NIFR are at 
relative rest. A notion of relative acceleration can be deter- 
mined in a covariant way [0]. This value is equal to zero. 
However, eqs.(||), strictly speaking, are not satisfied. 

If space and time are relative in the BLMP sense, then 
for observers in the IFR and NIFR the motion of the point 
masses forming the reference body, which are kinemati- 
cally equivalent, must be dynamically equivalent too (both 
in the nonrelativistic and relativistic sense). That is, if from 
the viewpoint of the observer in the IFR, the point masses 
forming the NIFR RB are at rest ( are not subject to the in- 
fluence of forces ),then from the viewpoint of the observer 
in the NIFR the point masses forming the RB of his frame 
are at rest too (are not subject to the influence of forces ei- 
ther ). In other words, if for the observer in the IFR the 
world lines of the IFR RB points are, according to eq. ([I]), 
the geodesic lines, then for the observer in the NIFR the 
world lines of the NIFR RB points are also the geodesic 
lines in his space-time, which can be expressed by eq. 
The differential equations of these world lines at the same 
time are the Lagrange equations of motion of the NIFR RB 
points. The Lagrange equations, describing the motion of 
the identical RB point masses in the IFR, can be obtained 
from the Lagrange action S by the principle of least action. 
Therefore, the equations of the geodesic lines can be ob- 
tained from the differential metric form ds — k dS, where 
k is the constant, dS = L(x,x)dt and L is the Lagrange 
function. The constant k = — (m p c) , as it follows from 
the analysis of the case when the frame of reference is iner- 
tial. It is equal to — (m p c) _1 . 

Thus, if we proceed from relativity of space and time in 
the BLMP sense, then the differential metric form of space- 
time in the NIFR can be expected to have the following 
form 

ds = — (m p c) _1 dS(x, dx). (3) 

In this equation S is the Lagrangian action describing (in an 
IFR) the motion of the identical point masses m p , forming 
the NIFR reference body. 

So, the properties of space-time in the NIFR are entirely 
determined by the properties of the used frame in accor- 
dance with the idea of relativity of space and time in the 
BLMP sense. 

Consider two examples of the NIFR. 

l.The reference body is formed by noninteracting elec- 
tric charges moving in a constant homogeneous electric field 
E. The motion of the charges is described in Cartesian co- 
ordinates by the Lagrangian 

L = -m p c 2 (l-V 2 /c 2 ) 1/2 +E ex, (4) 



where V is the speed of the charges. 

According to eq. (^) the space - time metric differential 
form in this frame is given by 

ds = dr\ — (wx/c 2 )dx° , (5) 

where dr) = (c 2 dt 2 — dx 2 — dy 2 — dz 2 ) 1 / 2 , is the metric 
differential form of the pseudo - Euclidean space - time in 
the IFR and w = eE/m is the acceleration of the charges. 

2. The reference body consists of noninteracting elec- 
tric charges in a constant homogeneous magnetic field H 
directed along the axis z. The Lagrangian describing the 
motion of the particles can be written as follows [|[ 

L = -m p c 2 {\ - V 2 /c 2 ) 1/2 - (m p n /2)(xy - xy), (6) 

where x = dx/dt, y = dx/dt and S7q = eH/2mc. 

The points of such a system rotate in the plane xy around 
the axis z with the angular frequency 

Lu = n Q [l + {Q Q r/c) 2 ]- 1 / 2 (7) 

where r = (x 2 + y 2 ) 1 ^ 2 . The linear velocities of the BR 
points tend to c when r — > oo. 
For the given NIFR 

ds = dt]+ (fi /(2c) (ydx- xdy). (8) 

In the above NIFR ds is of the form 

ds = T{x 1 dx), (9) 
where Tix, dx) = dr] + f a (x)dx a , f a is a vector field and 

drj = [r la(3 (x)dx a dx 13 } 1 / 2 

is the differential metric form of pseudo-Euclidean space- 
time of the IFR in the used coordinate system. T is a ho- 
mogeneous function of the first degree in dx a . Therefore, 
generally speaking, the space-time in NIFR is Finslerian iQ] 
with the sign-indefinite differential metric form 

4 Space and Time Measurements in 
NIFR 

For the 3 + 1 decomposition of space - time in noniner- 
tial frames of reference to 3-space and time we proceed 
from a covariant method which goes back to Ulman, Ko- 
mar, Dehnen and other authors [||. An ideal clock is a lo- 
cal periodic process measuring the length of its own world 
line 7 to a certain scale. For an observer in the NIFR the 
direction of time in the point x a is given by the vector of 
the 4-velocity £ Q of the BR point. 

The physical 3-space in each point is orthogonal to the 
vector The arbitrary vector £ Q in the point x a can be 
represented as follows: 

r = r + /?c, do) 



where £ are the spatial components and (3 is the function 
ofx a . 

Suppose any space vector £ in the point x a is orthog- 
onal to the vector ( a in the sense of the Finslerian metric 



(ii) 



where £ a is the co variant components of the vector £ Q , 
which are given by 



(12) 



Since T{x, ( a ) = 1 this vector is of the form 

C« = + fa = V a + f a , (13) 

where v a — dx a jdr\ is the 4-velocity of the reference body 
point in the IFR. By multiplying eq. ( |lC)| ) to the vector 
we find that f3 = ( a £ a and 



(14) 



where H % — ( a (p — Sp and 8^ is the Kroneker delta. 
Eq. ( jlOp yields for the vector £ Q = dx a : 

dx a = dx~ a + cd T ( a , (15) 

where dx° is the spatial components of the vector dx a and 

dr = c-\ a dx a (16) 

is the time element between the events in the points x a and 
x a + dx a in the NIFR. 

The metric form (^) and the spatial projection of the 
vector dx a lead to the following covariant form of the spa- 
tial element in the NIFR 



dL = {—rjafjdx dx )*/ 2 + f a dx 



(17) 



This covariant equation is the simplest and clearest in 
the coordinates system in which 



^=<5 7(%o) 1/2 



Indeed, in eq. ( |l7| ) 



-rjapdx dx^ — H af 3dx a dx 13 1 , 



(18) 



(19) 



where 



V^eCCaCp + CaC/3 + CfiC a , 

and ( Q = -riapC 13 

In the used coordinates system 

C" = v a dT)jds = \5%, 



(20) 



(21) 



where A = l/(^o + /o) an d = (^oo) 1 ^ 2 - The zero- 
component of the tensor H a p is 

= goo - 2r] oH' y o + fo) + {dr]/ds)(v Q + fof. (22) 

Since 

d V /ds = (1 + Uu 01 )- 1 = (7 ?00 ) 1/2 A (23) 

and A(^o + fo) = 1, the value of Hoo is equal to zero 
identically. 

The components 

Hik = r)ik + 2\rj t o{vk + fk) + 

(dn/da) 2 (yi + fi)(u k + f k ), (24) 

where v t = i]oi/{Voo) 1/2 - The spatial tensor H ik = ij ik 
with accuracy up to V/c, where V is the linear speed of the 
reference body points. 
We have also 



where 



f a dx = f a dx a , 



fa — faHp — fa - fo^{v a + fa)- 



(25) 



The zero-component of the vector f a is equal to zero iden- 
tically and the spatial-components are equal to /$ with ac- 
curacy up to V I c. 

Thus, in the used coordinates system the spatial element 
in the NIFR with accuracy up to V/ c is of the form 

dL = (~r lik dx l dx k ) 1/2 + fidx 1 = dl(l + (26) 

where dl — {—r\ ik dx l dx k ) x l' 1 is the Euclidean spatial ele- 
ment and k % = dx 1 /dl is a unit vector of the direction with 
respect to dl. 

The phase shift in the interference of two coherent light 
beams on a rotating frame was observed by Sagnac [gj. For 
a relativistic explanation of the effect it is usually postulated 
,that space-time in any frames of reference is pseudo - Eu- 
clidean [[Tot, [^J. The motion in NIFR is considered as the 
relative one in absolute pseudo- Euclidean space-time. 

However, for an isolated observer (in a "black box") in 
the rotating frame, who proceeds from the notion of space 
and time relativity in the BLMP sense, the observed aniso- 
tropy in the time of light propagation (which from his view- 
point contradicts the experiments of Michelson - Morley 
type) is not a trivial effect. It must have some "internal" 
physical explanation. 

A rigid disk rotating in the plane xy with angular ve- 
locity £1 around the axis z is approximately identical to the 
NIFR described in example 2 of Sec. 3. It follows from eq. 
( |26| ) that the spatial element in the NIFR is anisotropic. We 
will show that the speed of light in the noninertial frame of 
reference is anisotropic too. 

A triple of space basis vectors, necessary to compare 
the direction of a given vector from viewpoints of the NIFR 
and IFR, are not defined above in each point of the NIFR. 



A 



However, it does not prevent us from comparing the lengths 
of the vectors. In particular, we can find a dependency be- 
tween the speeds of a particle in the NIFR and IFR. 
The speed of the motion of a particle in the NIFR is 



V' =(-7 la(3 U a U^ 2 +f a U Q , 

where u a — dx a /dr is 4-velocity of the particle. 
The term under the square root is given by 



(27) 



= -r) a f3U a u - c 2 {dr]/ds) 2 + 2c 2 (rfr /dr), (28) 
where we have used the following equalities: 

C a dx a /dr = c, (29) 

( a = v a dr//ds, (30) 

( a dx a /dT = c (31) 

r = c~ 1 u^dx fl (32) 

The first term in eq. ( |28| ) for a photon is equal to zero and 
we obtain with accuracy up to V/c that {—rj a pu a u^) 1 / 2 = 
c. 

In the same approximation f a u a = / Q " a = cfik 1 . 
Thus, in the used coordinates system the speed of the 
photon in the NIFR with accuracy up to V/c is 



v' ph = c{l + f i k i ). 



(33) 



Consider a disk rotating with the constant angular ve- 
locity £1 around the z axis. Let r and 8 be the coordinates, 
defined by the equations 

x = r cos(tp) ,y = r sin(<p) , ip = 8 + fit . (34) 

In the coordinate system (r, 8, z, t) the space - time metrical 
differential form ds in the rotating frame is of the form 

ds = dr)- [flr 2 /(2c)]d8 - [n 2 r 2 / (2c 2 )]dx° . (35) 

where dr] is the pseudo - Euclidean metric form: 

d V 2 = [1 - (nr) 2 /c 2 ]{dx ) 2 - (dr) 2 - r 2 (d8) 2 



-2(r 2 Q/c)d8dx° - dz 2 



(36) 



In this coordinates system eq. ( |T8| ) is satisfied with ac- 
curacy at least up to V/ c. In virtue of equations ( 26 ) and 
( |33| ) the time of the motion of light from the point x l to 
x l + dx l is dL/vp = c~ 1 dl(l + 2fik l ). It follows from eq. 
( |36| ) that in the used coordinates system 

Jl 2c 2tt' 

For this reason the difference in the time interval between 
light propagation around the rotating disk in a clockwise 
and counterclockwise direction is Anr 2 Vt/c 2 , which gives 
the Sagnac phase shift [|[. Thus, the Sagnac effect for 
the isolated observer in the rotating frame can be treated as 
caused by the Finslerian metric of space-time in noninertial 
frames of reference. 



5 Inertial Forces 

Let us show that the existence of the inertial forces in NIFR 
can be interpreted as the exhibition of the Finslerian con- 
nection of space-time in such frames. 

According to our initial assumption in Section 3,the dif- 
ferential equations of motion in an IFR of the point masses, 
forming the reference body of the NIFR, are the geodesic 
lines of space-time in NIFR. These equations can be found 
from the variational principle 6 J ds = 0. The equations 
are of the form 



dC/ds + G a (x,C) = 0, 



(37) 



where ( a is the 4-velocity of the point mass, the world line 
of which is x a = x a (s), and 

G a (x, = T a pi CC + + CP dip' 1 ) /ds, (38) 

where 

p = dr)/ds = ( V apCC P ) 1/2 , B<$ = v a5 B Sf3 

and 

Bsp = dfp/dx 3 - dfs/dx?. 

In the Finslerian space-time a number of connections 
can be defined according to eq. ([?7|) [Q. In particular, this 
equation can be interpreted in the sense that in the NIFR 
space-time the absolute derivative of a vector field £ a (x) 
along the world line x a = x a (s) is of the form 



DC/ds = dC/ds + G%(x, dx/ds)^ 



(39) 



where 



G^{x,dx/ds) = Bf 3 +T^ 1 dx 1 /ds + pdp- 1 /ds 



Equations (39 ) define a connection of Laugvitz type [|7|] 
in space-time of the NIFR, which is nonlinear relative to 
dx a . The change in the vector due to an infinitesimal 
parallel transport is 



-G%(x,dx)^, 



(40) 



Consider the motion of a particle of the mass m p , in a 
NIFR, unaffected by forces of any kind in the laboratory 
(inertial) frame of reference. The differential equations of 
motion of such a particle can be found from the variational 
principle 5 J dr\ — 0. Since ds = dr] — f a dx a , the equa- 
tions of motion are 



Du a /ds = BZvP. 



(41) 



As an example, consider the nonrelativistic disk rotating 
in the xy plane about the z axis with the angular velocity £1 
The equations of motion ([37]) are 



dC/dt + Q, x r = 0, 



(42) 



where r — {x, y, z} and the coordinates origin coincides 
with the disk center. The absolute derivative (|9|) of a vector 
£ is given by 



D£/dt = d£/dt-Q x £. 



(43) 



and the equations of motion ( |4l| ) of the considered particle 
in the NIFR are 



Dv/dt = — fi x v, 

where v = {x y, z}. 

Next, for the 4-velocity u a we have 

u a = u a +C 



(44) 



(45) 



where u a is the spatial velocity of the particle in the NIFR. 
In the nonrelativistic limit eq. (45) can be written in the 
form 

v = v + (, (46) 

where v is the "relative" velocity of the particle and £ = 
V is the velocity of the disk point in the laboratory frame. 
Substituting (^) in Q), we find that 



Dv/dt = -D(/dt -nxv-QxC- 



(47) 



The value Dv/dt is an acceleration of the considered 
particle in the used NIFR found with the help of measuring 
instruments. The velocities field £ of the disk points is given 
by 

( = fixr. (48) 
Hence, along the particle path we have d£ /dt — ft x v 

and 

D(/dt = d(/dt-Q,x( = nxl!. (49) 
Thus, finally, we find from (Q) 

m p Dv/dt = —2m p (Q x v) ~ m p U x (fi x r). (50) 

We arrived at the nonrelativistic equations of motion of a 
point in a rotating frame jl2|]. The right-hand side of eq. 
( |5C| ) is the ordinary expression for the Coriolis forces and 
the centrifugal force in the rotating frames. See also [ |l3[ ] )■ 
Thus, in the nonrelativistic limit the Finslerian space- 
time in NIFR manifests itself in the structure of vector deriva- 
tives with respect to the time t. It should be noted that eq. 
d43| ) is considered sometimes in classical dynamics nomi- 
nally [ |l2| | just for the derivation of the inertial forces in the 
NIFRs. 

6 Relativity of Inertia 

A clock, which is in a NIFR at rest, is unaffected by acceler- 
ation in space - time of the frame. The change in rate of the 
ideal clock is a real consequence of the difference between 
the space - time metrics in the IFR and NIFR. It is given 
by the factor a — ds/dr\ from the equation ds = adrj. For 
the rotating disk of the radius R a = 1 — uj 2 R 2 /2c 2 which 



gives rise to the observed red shift in the well known Pound 
- Rebka - Snider experiments. 

We consider here another experimentally verifiable con- 
sequence of the above theory. 

Let p a = iripC dx a /dr] be 4-momentum of a particle 
in the IFR. Using 3 + 1 decomposition of space-time in the 
NIFR we have 

EC- (51) 



P a =p a 



From the viewpoint of an observer in the NIFR the spatial 
projection p a should be identified with the momentum, and 
the quantity cE with the energy £ of the particle. It is obvi- 

ous that E = C, a p a . 

Therefore, the energy of the particle in the NIFR is 



£ = m p Qc 2 £ a u a , 



(52) 



where Q — ds/di] — F(x, dx/drj). For the particle at rest 
in the NIFR u a = £" and we obtain 



£ = nipQc 



(53) 



Thus, the inertial mass m' p of the particle in the NIFR is 
given by 

m' p = Qm p (54) 

The quantity m' p coincides with the proportionality factor 
between the momentum and the velocity of the nonrelativis- 
tic particle in the NIFR. 

Since Q is the function of x a , the inertial mass in the 
NIFR is not a constant. For example, on the rotating disk 
we have 

m' p = m p /(l - n 2 r 2 /2c 2 ), (55) 

where ft is the rotation angular velocity and r is the distance 
of the body from the disk center. 

The difference between the inertial mass m p q of a body 
on the Earth equator and the mass m v p ° l of the same body 
on the pole is given by 



{m e p q - wP p ol )/ml° l = 1.2 • 1(T 12 



(56) 



The dependence of the inertial mass of particles on the Earth 
longitude can be observed by the Mbssbauer effect. Indeed, 
the change AA in the wave length A at the Compton scatter- 
ing on particles of the masses m p is proportional to rn p 1 . If 
this value is measured for gamma - quantums with the help 
of the Mossbauer effect at a fixed scattering angle, then af- 
ter transporting the measuring device from the longitude (pi 
to the longitude if2 we have 



(AA)" 1 - (AA) 



•p'l 



(AA)^ 1 
where 9 = 1.2 ■ 1CT 12 . 



Q[cos 2 ( l p 1 )-cos 2 (^ 2 )}, (57) 



7 Gravitation in Inertial and Proper 
Reference Frames 

Consider a frame of reference whose reference body is formed 
by identical material points m p moving under the effect of 
the field tp a p. These frames will be called the proper frames 
of reference (PRF) of the given field. Any observer, located 
in the PFR at rest, moves in space-time of this frame along 
the geodesic line of his space-time. This implies that the 
space-time metric differential form in the NIFR is given by 
eq. (|j) where S is the action describing describing the mo- 
tion of particles forming the reference body of the NIFR. 

Now suppose jl7[ | that in pseudo- Euclidean space-time 
gravitation can be described as a tensor field ip a p( x ), and 
the Lagrangian describing motion of a test particle with the 
mass m p is of the form 

. a . (3 , 

L = -m p c[g a p(ip) x x } ' , (58) 

where x = dx a jdt and g a p is the symmetric tensor whose 
components are the function of ip a p . 

According to (]3|) the space-time metric differential form 
in the PFR is given by 

ds 2 = g a p{tp) dx a dx 13 (59) 

Thus, the space-time in the PFR is a Riemannian with 
the curvature other than zero. Viewed by an observer in the 
IFR, the motion of the test particle forming the reference 
body of the PFR is affected by the force field ip a p. But 
the observer located in the PRF will not observe the force 
properties of the field ip a p since he moves in space-time of 
the PRF along the geodesic line. For him the presence of 
the field ip a p will be displayed in another way — as space- 
time curvature differing from zero in these frames, e.g. as 
a deviation of the world lines of the neighbouring points of 
the reference body. 

For example, when studying the Earth's gravity, a frame 
of reference fixed to the Earth can be considered as an iner- 
tial frame if the forces of inertia are ignored. An observer 
located in this frame can consider motion of the particles 
forming the PRF reference body in flat space-time on the 
basis of eq. ( |58| ) without running into contradiction with ex- 
periments. However, the observer in the PFR (in a comov- 
ing frame for free falling particles) does not find the Earth 
gravity as some force field. If he proceeds from the relativ- 
ity of space-time, he believes that point particles, forming 
the reference body of his reference frame, are the point of 
his physical space. They do not affect a force field and their 
accelerations in his space-time are equal to zero. In spite of 
that, he observes a change in the relative distances of these 
particles. Such an experimental fact has apparently the only 
explanation as non-relativistic display of the deviations of 
the geodesic lines caused by space-time curvature. So, we 
observe an important fact that only in proper frames of ref- 
erence we have an evidence for gravitation identification 
with space-time curvature. 



Thus, we arrive at the following hypothesis. In iner- 
tial frames of reference, where space-time is pseudo - Eu- 
clidean, gravitation is a field ip a p- In the proper frames 
of reference of the field ip a p, where space-time is Rieman- 
nian, gravitation manifests itself as curvature of space-time 
and must be described completely by the geometrical prop- 
erties of the letter. 

If this possibility really takes place in nature, then it will 
remove an isolation of the geometrical gravitational theo- 
ries from the theories of other fields. 

Of course, eq. (^) refers to any classical field. For in- 
stance, space-time in the PRF of an electromagnetic field is 
Finslerian. However, since ds depends on the mass m p and 
charge e of the point masses forming the reference body, 
this fact is not of great significance. 

It should be noted that the geometrical theory of gravi- 
tation in the PFR is not identical to Einstein's theory. Gravi- 
tational equations should be some kind of differential equa- 
tions for the function ip a p or g a p (-0), which are invariant 
under a certain set of gauge transformations of the poten- 
tials ip a p. Since g a p = g a p (-0), the Einstein equations are 
the equations both for g a p and for ip a p. Under the transfor- 
mation i\) a p — > ip a g the quantities g a p (ip) undergo some 
transformations too, and, as a consequence, the equations 
of the test particle motion resulting from eq. (|58|) and the 
Einstein's equations do not remain invariant. 

The equations of motion resulting from eq. (|8|) are at 
the same time the equations of a geodesic line of the Rie- 
mannian space-time V n of the dimensionality n with the 
metric tensor g a p (iji). That is why if the given gauge trans- 
formation ip a p — > ip a g leaves the equations of motion in- 
variant, then the corresponding transformation g a p — > ~g a p 
is a mapping V — ► V of the Riemannian spaces leaving 
geodesic lines invariant, i.e. it is a geodesic, (projective) 
mapping. Let us assume that not only eq. ( Jssj ) but also the 
field equations contain ip a p only in the form g a p (ip), then 
it becomes clear that the gauge-invariance of the equations 
of motion will be ensured if the field equations are invari- 
ant with respect to geodesic mappings of the Riemannian 
space V n . Thus, if we start from eq. (|58|), then the gravita- 
tional field equations as well as the physical field character- 
istics must be invariant with respect to geodesic (projective) 
mappings of the Riemannian space-time V n with the metric 
tensor g a p (ip). 

Simplest equations of that kind are analyzed in the next 
Section. 

8 The geodesic-invariant equations of 
gravitation 

In accordance with the basic principles the space-time in 
the PRF of the gravitational field ip a p is the Riemannian 
V n of dimension n=4 whose metric tensor is defined up 
to geodesic (projective) mappings g a p — > g a g, generated 
by the gauge transformations tp a p — > ip a p- The geodesic 



transformations of the metric tensor g a p are given by Levi- 
Chivita equations [|l9[, 



where 



where 



d 



(60) 



1 



•IP 



V a g pi = dg Bj /dx a - r? * 



2(n + 1) dx° 
.9 = det Iso^l 

The Christoffel symbols and the curvature tensor also 
do not remain invariant. In particular, the Christoffel sym- 



(61) 



bols are transformed as follows [BOJ: 



where <fi a (x) is a vector field. 

However, some objects, which are invariant under geodesic 
(projective) mappings of the space V n , also can be defined. 
Just these gauge-invariant objects have a physical sense in 
the theory under consideration. 

The simplest gauge-invariant object is the Thomas sym- 
bols [01: 



TT 7 

LL a0 



U f] L Ctfi 



(62) 



The Thomas symbols are transformed as a tensor only 
with respect to the projective coordinate transformations. 
However, by replacing the ordinary derivatives to the co- 
variant ones for the metric da 2 , a tensor B^g can be re- 
ceived. This object also can be written as follows 



BZp = Kp-n aP ,, 

o 7 

where H a3 are Thomas symbols in the E n : 



(63) 



o7 o7 

n Q/3 =r Q/3 -(n + i) 



SZ r e/3 +s} r ea 



(64) 



This geodesic - invariant tensor will be named the strength 
tensor of a gravitational field. Note that the equality £? 7 7 = 
is satisfied identically. 



According to eq. (58 ) the differential equations of mo- 



tion of the test particle are given by 

d 2 x^ 



ds 2 



dx^dx^_ 
^ ds ds - 



(65) 



where T 7 ^ are the Christoffel symbols in V n . Following 
[ p2] ] we will show how one can define a geodesic - invariant 
connection and curvature in the spaces under consideration. 

Let us define a scalar parameter p on the geodesic lines, 
that remains unaltered by geodesic mapping of V n , by means 
of a differential equation 



{p,s} = -2eT° a 



dx a dx 
ds ds 



(66) 
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\ ds 



(67) 



e is a nonzero constant and T" g are a given function of x a 
symmetric in low indices: 



p0 _ -pO 
1 a/3 — 1 0a- 



(68) 



By eq. (f56|) the parameter p is defined as the function of s 
up to linear fractional transformations. 

Since the parameter p must be a scalar, the object T° g 
is the components of the covariant tensor of rank 2. 

Let T a p and ds be the components of T° g and the line 
element ds, respectively, after some geodesic mapping of 
space-time V n . Then the new geodesic equations are given 
by 

d 2 x^ - 7 dx a dx? n 

^F +r ^r^r = °- (69) 

On the other hand, after a geodesic mapping eqs.(|65|) are 

d 2 x^ s" dx a 7 dx a dx 13 _ 
~W + (s') 2 ds + ds ds ' 

where s' = ds/ds , s" = d 2 ~s/ds 2 , and these equations 
must be identical to d69h. With regards to (plh, this yields 



ds/ds = exp ^—2 J (f> a dx a 

Then, by putting 

r _-. dx a dx 13 
{p,s\ = -2eT a8 = = 
ds ds 

and by taking notice of 

{p, s} = (ds/ds) 2 [{p, s} - {§, s}} , 



(70) 



(71) 



(72) 



we find (since p must be the invariant under geodesic map- 
pings) the equations of transformation of J?° B under geodesic 
mappings: 



r° -r° , 



-i r i / ^7 \ 

2 ^da; 7 + dxP ) 



(73) 



Now on every geodesic we define a gauge variable ("5 
th coordinate") by substitution 

A 1 ds 
2^ = --log-. (74) 
2e dp 

Let s be another parameter on the geodesic and 

3f* = _llog* (75) 
2e dp 

is the gauge variable corresponding to the transition p — > s. 
Then it follows from ( |70| ) that for any path 



x° = x° + e" 1 / (f> a dx a + C, 



(76) 



o 



where integration is performed from an arbitrary fixed point 
q on the geodesic along the curve and C is the arbitrary 
constant. 

Using the relation 



{v,u} 



we have 



{s,p} = -2eT° a0 



{u,v}, 



dx a dx@ 



(77) 



(78) 



dp dp 

On replacing in this equation and the geodesic equations 
(M) the derivations ds /dp, d 2 s/ dp 2 ,d 3 s/ dp 3 by x° , dx° / dp, 
(Fx /dp 2 according to the definition of the gauge variable 
( J74| ) we find that the geodesic equation ( p5| ) and the equa- 
tion ( f78| ) which yields the definition of the parameter p can 
be written as 



dp 2 
d 2 x 



■fa' 



dx@ dx 1 
dp dp 



^ dx A dx a q 
dp dp 



(79) 



2,4 /dj ,4 

e 



dp 2 \ dp 

respectively. These equations can be united as an equation 
of a geodesic in 5-dimensional space-time 



d-K 2 



BC 



dx B dx c 
d-K dn 







(81) 



The capital indices run from to 4 and r% = T$ B 
M5g where 8g is the Kronecker symbols. 

Suppose the functions <fi a satisfy the conditions: 



dfa/dx? - d<t> p /dx a = 0. 
If we consider the transformations 



and 



x° = x a 



,dx a + C, 



(82) 



(83) 



(84) 



where <fi a dx a is the exact differential of a function of x a , 
as admissible coordinate transformations in the n + 1 di- 
mensional manifold M.5, then eqs. (|Tj) can be regarded as 
the differential equations of the geodesic lines in homoge- 
neous coordinates of projective geometry. A theory of the 
projective connection in such a way has been considered by 
@] and other authors in 1921-1937. 
The object components Tg C are transformed as fol- 
lows: 



is transformed as the tensor relative to the coordinate trans- 
formations in M§. The components of the tensor TZg CD 
vanish identically save ftig™ and 



/C /37<5 



Dtt 

n f3j5 " 



where R^ s 



0S u y (87) 
is the curvature tensor of the affine connection 



in V4. It has the following properties 





V a 


= 0, 


(88) 






= 0, 


(89) 






= 0. 


(90) 



It follows from eq. (|37|) that the contracted tensor is 
given by 

ft/37 = -R/37 + e(n - 



(91) 



It does not change with geodesic mappings of V4. 
The equations 

n M = (92) 

are the simplest geodesic - invariant generalization of the 
Einstein vacuum equations. 

Depending on a choice of the object and e we ob- 
tain a specific variant of the theory. In this paper we will 
assume that e = 1 and 



1 a/3 



1 



(n + 1) 
[dQ a /dx + dQ /dx a 



(93) 



where Q a = r* 3 



/3a 



o/3 

r 



/3a 



The object defined in this way 



has the required properties under geodesic mappings (73) 
Equation (p2|) can be written in the form 



pa 







(94) 



where the semi-colon denotes a covariant derivative in the 
Pseudo - Euclidean space-time E4. These equations were 
proposed first in [ 18]] from another viewpoint. 

Equations (|94 ) are the system of the differential equa- 
tions for the geodesic- invariant tensor Bp (or for the func- 
tions g a p) which are defined up to arbitrary geodesic map- 
pings. The coordinate system is defined by the used mea- 
surement instruments and is given. The equations do not 
contain the functions tp a p explicitly. 

The simplest way of obtaining equations for ip a p is to 

set 



r 



BC 



^ F dx E dx D 
ED dx B dx c 



d 2 x F \ dx J 
dx B dx c ) dx 1 



(85) 



The components of Tg C are coefficients of the projec- 
tive connection [E2|. The object 



T^bcd = dn c /dx D - dT* D /dx 



'^MD^BC ^MC^BD (86) 



TDOt 

B 0i 



V^-Cn + l)" 1 ^^ 



(95) 

where V Q is the covariant derivative in flat space-time. The 
identity £?2 7 = is satisfied as it is to be expected accord- 
ing to the definition of the tensor B^. 

Then, at the gauge condition V <T ?A ( j 7 = eq. ( p4j ) are 
given by 



(96) 



where □ is the covariant Dalamber operator in pseudo - Eu- 
clidean space-time. 

It is natural to suppose that with the presence of matter 
these equations are given by 



(97) 



where k = 87rG/c 4 , t a p = x~ 1 \7 a ijj ai 'V y vp a f3 and T a p 
is the matter tensor of the energy-momentum. 
Obviously, the equality 



(98) 



is valid. Therefore, the magnitude t a p can be interpreted as 
the energy-momentum tensor of a gravitational field. [] 



Spherically-Symmetric 
tional Field. 



Gravita- 



Let us find the spherically symmetric solution of eqs.(^4j) in 
an inertial frame of reference, where space-time is supposed 
to be pseudo - Euclidean. 

Because of the gauge (geodesic) invariance, additional 
conditions can be imposed on the tensor g a p. In particular 
0181, under the conditions 







(99) 



eqs. (|94|) will be reduced to the Einstein vacuum equations 
Ra/3 = 0, where R a p is the Ricci tensor. Let us choose a 
spherical coordinate system. Then, if the test particle La- 
grangian is invariant under the mapping t — > —t, the funda- 
mental metric form of space-time V4 can be written as 



ds 2 



-A(dr) 2 - B[{d6) 2 + sm 2 (6)(dtp) 2 



-C(dx 



0\2 



(100) 



where A, B and C are the functions of the radial coordi- 
nate r. Proceeding from the above- stated, we shall find 
the functions A, B and C as the solution of the equations 
system 

R a p = (101) 



and 



Q a = 0, 



(102) 



which satisfy the conditions: 



lim A = l, lim (B/r 2 ) = 1, lim C = 1. (103) 



1 It should be noted that, when we introduce it in some way, we can- 
not be sure apriori that the equation for ij) a p yields all solutions of the 
equations for . We may introduce a potential i/) a( 3 also in another 
way. 



The equations R±± = and i?oo = can be written 

as 

BL 2 - 2CL X = (104) 



and 



B'C - 2BL 2 = 0, 



(105) 



where 

Li = #1212 = B"/2 - (B') 2 /m - B'A'/{AA) , 

L 2 = Rwio = -C"/2 + C'(AC)'/(4AC) 

and the differentiation of A, B and C with respect to r is 
denoted by the primes. 



The nonzero eqs.( |102[ ) yield 



B AC = 



(106) 



First, the combination of eqs. (|104j) and ( |105| ) yields 

ACLt - B'C = 0, 



i.e. 



2B" - {B'f/B - B'(AC)'/(AC) = 0. 
Also, taking the logarithm of eq. ( |106| ) we obtain 

(AC)' /(AC) = 4/r - 2B'/B = 0. 



(107) 



(108) 



Equation (107) then becomes 

2B" + (B'f/B - AB'/r = , 



or 



v! /u = 4/r, 



where u = (B') 2 B. By using (103) we find 

B = (r 3 +JC 3 ) 2 / 3 , 
where JC is a constant. 



(109) 



(HO) 



Next, from eq. (105) we find by using eqs. (103), 

C=l-Q/B 1 ' 2 , (111) 
where Q is a constant. 



Finally, we can find the function A from eq. ( 106 ). 
Thus, in the spherical symmetric coordinate system the 
following functions A, B and C are obtained: 

A = (f') 2 (l - Q/f)-\ B = f 2 ,C = l- Q/f, (112) 

where 



/ = (r 3 +/C 3 ) 1 /3 



The nonzero components of the tensor B^p are given 



and /' = df/dr. 
by 



1 r\ 



rr 2 A ' 



^ " 2 rB 



1 2rA - 5' 

2 I ' 

1 2rA - B' 2 

= — t: -, sin t 

2 A 

1C „, 1 C" 



2 A ' 

r><t> _ 



2 C ' 
2B - B'r 

Vb 



(113) 



If K = 0, then / = r, and the solution ( |1 12[ ) coin- 
cides with the Droste-Weyl solution of the Einstein equa- 
tions which is known as the Schwarzshild one 
it coincides with the originally Schwarzshild solution |[Z6||. 
However, it is important to understand th at fro m the point 
of view of the considered theory, solution ( 112 ) is obtained 
in a given (spherical) coordinate system, defined in pseudo 
- Euclidean space-time, and that different values of the con- 
stants Q and fC yield different solutions of equation ( p4| ) in 
the same coordinate system. 

The equations of the motion of a test particle resulting 
from Lagrangian (p8b is given by 



'x a + (T? jl -c- 1 r ( / jj x a )x x~< = 0. (114) 

In the nonrelativistic limit x = — c 2 Fq , where Tq = 
C/2A = r^C'IfC. Therefore, to obtain the Newton 
gravity law it should be supposed that at large r the func- 
tion / w r and Q — r g — 2GM/c 2 is the classical 
Schwarzshild radius. 

At a given constant Q allowable solutions are obtained 
by change of the arbitrary constant JC. In particular, if we 
setting K, = the fundamental form of space-time V4 co- 
incides with the Droste-Weyl solution of the Einstein equa- 
tion [ |28[ ] (it is commonly named the Schwarzshild solution) 
which have an events horizon at r = r g : 



ds 2 = - 



+ (!-!£) dx° 



(115) 



If we setting K = Q, the solution coincides with the orig- 
inal Schwarzshild solution [ p6| ] which have no the event 
horizon and singularity in the center: 



ds 2 



f 2 dr 2 
+ (l- r JL)dx° 2 , 



f 2 [d0 2 +sm 2 0d(p 2 



(116) 



where / = ( 7^ + r 3 ) 1 / 3 . 

These solutions are related to the same coordinate sys- 
tem and are different solutions of the gravitation equations 
under consideration. 



2 A discussion of the difference in these solutions are given in | 



Of course, a reader can say: however, we can obtain the 
Droste-Weyl solution from the original Schwarzshild one 
by a coordinate transformation. Supposing it is so. (There 
is an alternative point of view [?]). However, in this case 
we must transform also space and time intervals to the new 
co-ordinates which, in contrast to the spherical ones, have 
no sense of values measured by rulers and clock. After 
that, of course, we obtain the same physical results as in 
the spherically co-ordinate. ( Like classical electrodynam- 
ics in arbitrary co-ordinates). Therefore, since at JC = Q 
the solution of our equation have no the singularity in the 
center and events horizon, it does not contain theirs and in 
others coordinate systems. 

We can argue that the constant Q = r g . Indeed, con- 
sider the 00-component of eq. (p7|). Let us set T a p — 
pc 2 u a up 1 where p is the matter density and u a is the 4- 
velocity of matter points. At the small macroscopic veloci- 
ties of the matter we can set uq = 1 and u% = 0. Therefore, 
the equation is of the form 



where \ 
tensor 



= X(pC 2 + t Q0 ) (117) 

4ttG/c 4 and t o is the 00-component of the 



(118) 



tap = X 1 BZa B 1 - 

Let us find the energy of a gravitational field of the point 
mass M as the following integral in the pseudo - Euclidean 
space - time 

£= ftoodV, (119) 



resulting from the above solution, where dV is the volume 
element. In the Newtonian theory this integral is divergent. 
In our case we have: 



too = 2x 1 BL B" 



c 4 Q 2 
00 " 01 ~ SirG K 



(120) 



and, therefore, in the spherical coordinates, we obtain 

r 1 q 2 c 4 



8 ttG 



where 



and 



B(z, w) 



00 t z-l 



dt 



(121) 



(122) 



(123) 



/o (l + t) 2+tu 
is B-function. Using the equality 

*M = ^4, (124) 
1 (z + w) 

where T is T-function we obtain J = 4tt/7C, and, therefore, 



Q Qc 4 Q 
£ = — — = — Mc 2 
K 2G K 



(125) 



We arrive at the conclusion that at JC ^ the energy of 
the point mass is finite and at JC = Q the rest energy of the 
point particle in full is caused by its gravitational field: 

£ = Mc 2 . 

The spacial components of the vector P a = to a are 
equal to zero. 

Due to these facts we assume in the present paper that 

K 



Q 



and consider solution (112) in the spherical 



coordinates system at the used gauge condition as a basis 
for the subsequent analysis. 

10 Orbits of Non-Zero Mass Parti- 
cles. 

The equations of motion of a test particle of a non-zero 
mass in the spherically symmetric field resulting from eqs. 
(M) are given by | 



(126) 



r = (c 2 C/A)[l~(C/E)(l + r 2 J / B% 



^ = c CJr g /(BE) (127) 

where (r, ip, 8) are the spherical coordinates (9 is supposed 
to be equal to n/2), r = dr/dt, ip = dip/dt , E = 
E/(mc 2 ), J = J/(amc), E is the particle energy, J is the 
angular momentum. Let u = 1//, where / = (1 + r 3 ) 1 / 3 
and r = r/r g . Then the d iffer entia l equation of the orbits, 
following from eqs. ( |l26| ) and (127) can be written as 



{dm/d^f = g{u) 



(128) 



where 



G(u) = u 3 - u 2 + u/ J" + (E~ - 1)/ J" 



eq. ( |l28[ ) differs from the orbit equations of general 
relativity j27|] by the function / instead of the function r. 
Therefore, the distinction in the orbits becomes apparent 
only at the distances r of the order of 1 or less than that. 

1 2 

Setting r — in eq. (|126j) we obtain E = 7V(r), where 



Af(r) 



a 



!//)(! + J 2 //) 



(129) 



is the effective potential ||2^]. Fig.([l]) shows the function 

TV = Af(r). 

The function jV(r) differs from the one in general rela- 
tivity in two respects: 

1. It is defined at every point of the interval (0, oo) 

2. It tends to zero when r — * 0. 

Possible orbit types can be shown by the horizontals 
E = Const. Two types of the orbits have peculiarity in 
comparison with the Einstein equations. The horizontals 
placed above the maximum of the curve Af(r) show the 
particles orbits which begin in the field center and end in 
the infinity. In other words, for each value of J there exists 




Figure 1 : The effective potential J\f — J\f(r) of the particles 
with J = 3, 3 1 / 3 , 0. (Curves 1, 2 and 3, correspondingly). 



such a value of E for which the gravitational field cannot 
keep particles escaping from the center. The events horizon 
is absent. Fig.|| shows an example of the orbits at J = 1.99 
andS = 1. 




o 



Figure 2: The orbit of the particle with J = 1.99 and E = 
1. 

The horizontals placed between Y-axis and the curve 
JV(r) can show particles orbit kept by the gravitational field 
near the field center. 

It follows from eqs. (126) and (127) that the velocity 



of a test particle freely falling to the point mass M tends to 
zero when r — > 0. The time of the motion of the particle 
from some distance r = tq to r = is infinitely large. 
We can say that the spherically symmetric solution has no 
physical singularity. 

The points of the minimum of the function J\f(r) show 
stable closed orbits, the points of the maximum show insta- 
ble ones. The minimum of the function 7V(r) exists only 
at J > \/3 which corresponds to the value of the func- 
tion /(r) > 3. Therefore, stable circular orbits exist only 
at r > r cr , where r cr = \/26 r g w 2.96 r g . The or- 
bital speed of the particle with r — r cr is equal to 0.4c. 
At r < r cr unstable circular orbits can exist. At J — > oo 



the location of the maximums tends to / = 3/2. There- 
fore, the minimum radius of the instable circular orbit is 
r m in — 1-33 r g . (In general relativity it is equal to 1.5 r g ). 
The speed of the motion of a particle on this orbit is equal 
to 0.51c. The binding energy E — 0.0572, just as it occurs 
in general relativity. 

The rotation frequency uj r = <p of the circular motion 
will be 

cu r = [cQ -l)/{fr g )](J/E) (130) 

In a circular motion r is the constant and, therefore, the 
function JV(r) has the minimum. Consequently, from the 
equation dj\f / dr = we find 



J 2 = / 2 /(2/-3) 



Using ( |126| ) we have at r = 



£T = 2(7-l)/[7(2/-l)] 



(131) 



(132) 



Equations (129)-=- d!31h yield 



u r = rV V(V2/ 3 / 2 ) 
. Hence, the circular orbits have the rotation period 

T = 2V27rc-V 3 (l + r 3 ) 1 / 2 (133) 

( 3 rd Kepler law). In comparison with general relativity the 
change in T is 2% at r = 3 and 20% at r = 1.33. 

Consider the apsidal motion. For ellipsoidal orbits the 
function Q(u) has 3 real roots u% < U2 < E 3 ^7p. The 
apsidal motion per one period is 



where 



A = 



Sip = 2|A| - 2tt, 

[G{u) 1/2 du. 



Consequently, [ 3 1 1 

A = 2(u 3 -ui) 1/2 W2,9), 

where 

q = [(e 2 -e 3 )/(e 3 - ui)} 1/2 



and 

T(ir/2,q) 



tt/2 



(1-q 2 sin 2 C3)) 1/2 d/3. 



(134) 
(135) 

(136) 
(137) 

(138) 



Let us introduce (by analogy with general relativity) the 
following notations: 

d = (1 - e)/p, E 2 = (1 + e)/p, E 3 = 1 - 2/JJ, (139) 

where the parameter p at r — > oo becomes the focal param- 
eter p divided by r g . 

At r > 1 the value of q ps (Te/p) 1 / 2 <C 1 and, there- 
fore, 

F(ir/2, q) = (tt/2)(1 + q 2 /A + 9g 2 /64 + ...). (140) 



Using eqs. ( 135 ) and ( 138) we find with accuracy up to 

1/p 2 

Aip = 3n/p + (7r/(8p 2 )(-e 2 + 16e + 54 + ...). (141) 

For the orbits of Mercury or a binary pulsar (such as 
PSR 1913 + 16 ) the value of E differs very little from the 
value of r g /r. Consequently, the values of p = 2/(Ei+E2) 
and e = (7*2 — Ei/(E2 + El) differ very little from the val- 
ues of pj r g and the e. Hence, their apsidal motion differs 
very little from the general relativity prediction. Even, for 
example, at p — 10 and e = 0.5, the difference in Aip is 
about 6 • lO^^rad/year. 



11 Photon Orbits. 

The equations of motion of a photon in the spherical sym- 
metric field are given by | 



r = (c 2 C/A)(l-C b 2 /B), V = cCb/B, (142) 

where b is an impact parameter. 

The differential equation of the orbits can be written as 



where 



(du/dip) 2 = Qi(u), 



0i(E) = u i ~u 2 + l/b 



(143) 



and b = b/ r g . 

Setting r = in eq. ( |143[ ) we obtain b = //(l - 
l//) 1 / 2 where / = f/r g . Fig || shows b as a function 
of r, i.e. the location of the orbits turning points. 
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Figure 3: The impact parameter b as the function of r 

The function b(r) is defined at all r > 0. The mini- 
mal value of b, i.e. 6 m j„, is equal to 2.6. It is reached at 
Tmin — 1.33. The motion of photons can be shown by the 
horizontal 6 = Const. There are two types of orbits which 
have the peculiarities due to the lack of the events horizon. 
The curves placed below the minimum of the curve b(r) 



show that the attracting mass cannot keep a photon escap- 
ing from the center at the parameter b < b m i n .The orbits 
of this type also show the gravitational capture of the pho- 
ton. The photon finishes at the field center, unlike general 
relativity, where it ends on the Schwarzshild sphere. 

The angle of the light deflection at the distances close 
to r = 1.33 is given by 



d 2 r A dx° dr 
ds 2 2A ds ds 

{dip 
2A I ds 



2A 

2 



dr 
ds 

C_ 

2A 



A dx° dr 
2A ds ds 

^ 2 
ds 



= 0, (146) 



ln(4.021/,5 2 ), 



(144) 



where S = 2(f mm - 3/2) and f mm = (1 + rt m ) 1/3 . 
It differs very little from the one in general relativity in a 
weak field. (See also [p[ ). 
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(147) 



12 Spherical-symmetric solution in 
general case 

Properties of the gravitational field, by definition, display 
themselves exclusively owing to its influence on the motion 
of test particles. Let us show that the motion of test parti- 
cles described by Lagrangian ( |58| ) is insensitive to the fact 
whether A, B, C are depending of tim e , or not. 

Suppose, in the metric form ( IOC ) the coefficients A, B, 
C are the functions of r and x°. The Christoffel symbols 
for this form are given by 



u 2C' 

rr 2C , 
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In these formulas we denote by prime the partial derivative 
with respect to r and by point - the partial derivative with 
respect to a; . Then the equations of the motion can be writ- 
ten in the form 



2^,0 



d 2 X 

ds 2 



A_ 

2C 



C_ 

2C 

dr 
ds 



dx° 
ds 



C dx° dr 
C ds ds 
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2C 
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0. 



(145) 



fp 
ds 2 



B dx° dp 
B~ds~~ds~ 



B' dr dp 
~B~d~s~ds~ 



= 0. 



(148) 



In these equations the prime and point denote the partial 
derivatives with respect to r and da; , respectively. 

Now let us find the motion integrals. 

Equation (147) has the solution 9 = tt/2. For this 



reason it can be assumed that the orbits are in the plane 



9 = tt/2. Equation (148) can be written in the form 



d 2 p 1 dB dp 
ds 2 B ds ds 



0. 



or 



where 



d 
ds 



ds 



ln(B) 



= 0, 



(149) 



dB _ dB dx a dB dr 
ds dx° ds dr ds 

is the total derivative of the function B (r, x°) with respect 
to s along the world line of the particle. It yields the first 
integral of the motion: 



B d f = J, 
ds 



(150) 



where J is a constant. 
Using the relation 



dr 



B 



dp 



rfs J \ds 
that follows from the identity 
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d6 



= 0, (151) 
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(152) 
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at 9 = tt/2, equation ( |145| ) can be transformed to the form 

d 2 x° 1 dC dx° 



0. 



(154) 



dt 2 C ds ds 
where dC/ds is the total derivative of the function 
C (r, a; ) with respect to s. That yields the second integral 
of the motion 



da; 
ds 



where A is a constant, eq. (146), after being multiplied by 



the factor 2 Adrj ds can be written in the form 



d_ 

ds 



A 



A 



dx° ( dr 



ds V ds 
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ds I ds \ ds 1 ds 



(156) 



and, after that, by using relation (151 ), in the form 



d 
ds 



A 



dB (dipY dC fdx° 



ds V ds 



+ 



ds 



= 0, 
(157) 

where dB/ds is the total derivative of the function 
B (r, a; ). Taking into account (150) and (155), we obtain 



d_ 

ds 



A[ dr 

ds 



ds\B ds\C 



which yields the third integral of the motion: 

, fdr\ 2 A 2 J 2 ^ 



= o, 

(158) 
(159) 



where E is a constant. 

Now it can be demonstrated that the same integrals of 
the motion can be obtained from Lagrangian (|58|) for the 
static field. At 6 = 7r/2 it can be written in the form 



2 2 

r A + ip B 
L 



L 



m 2 c 4 



But it follows from eq. (160) that 
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Ar 
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C 



= Const. 



L 



L 



(155) Consequently, C/L — Const. However, 



ds 



(163) 



(164) 



Thus, we arrive at equation (155). 

The se cond integ ral o f the motion can be fou nd from 
equations (162) and (164) and coincides with (15C). 

To obtain the third integral of the motion we start from 
the identity 



9af3- 



dx a dx 13 



c, 



(165) 



ds ds 

where ( or the particle with the mass rn p ^ and £ = 
for the particles with the mass m p = 0. Substituting (165) 
d(p/ds and dx°/ds from eqs. (163) and (164) into this 
equation, we obtain the equality which coincides with eq. 

Since the properties of a gravitation field are defined by 
their influence on the motion of the particles, a solution of 
any correct gravitation equations for the spherical symmet- 
ric field, based on Lagrangian ( |160| ), must be static. 

Using the above results, we can also find gravitational 
field inside a spherically - symmetric matter layer. In order 
to reach a coincidence of the motion equations of the test 
particles in the nononrelativistic limit with the Newtonian 



ones, the constant r g in eq. (Ill) in that case must be set 



equal to zero. Therefore, the spherically - symmetric matter 
layer does not create the gravitational field inside itself. 

13 Equilibrium Configurations with 
Large Masses 



L 



. 2 



m p eT c z C - Ar — Bip 



1/2 



(160) 



where A, B and C are the functions of r only. Since L 
does not depend on x° and <p, there are two integrals of the 
motion: 



■dL .dL _ _ M 
r — - + ip — : L = Const; 

dr dip 

— - = Const. 

dip 

The first of these equations is of the form 



(161) 
(162) 



In this Section we consider one of the most interesting 
consequences of gravitation equations ([w]) - the possibility 
of the existence of compact configurations of degenerated 
Fermi- gas with very large masses. 

The radial component of the gravity force affecting a 
test particle at rest in the spherically - symmetric field is 
given by | 



F' 



-mB, 



GrripM 



oo 



( r 3 + r 3)l/3 



(166) 



It follows from this figure that \F\ reaches its maximum 
at the distance r of the order of r g and tends to zero at r — > 
0. 



i c 




2 4 ,/r 6 8 10 



the nonrelativistic electrons: 
1.3 • 10 42 £/, .R = 2.3 ■ 10 12 cm, 

the relativistic electrons: p = 
10 40 g, R=1.Z- 10 n cm, 

the nonrelativistic neutrons: 
3.9 • 10 35 g,R= 1.6 • 10 7 cm. 



10' 



10 5 g/cm 3 , M 
g/cm 3 , M 
10 14 g/cm 3 M 



2.3 



The reason of the two types of configurations existence 
can be seen from fig. |[ where for p = 10 15 g/cm 3 , K = 
5 • 10 9 and T = 5/3 the plots of right-hand and left-hand 
sides of eq. (168) against mass M are shown. 



Figure 4: The modulus of the gravitational force (arbitrary 
units) as the function of the distance from the center (r/r g ). 



It would therefore be interesting to know what masses 
of the equilibrium configurations can exist if gravitational 
force is given by eq. ( 166 ). To answer this question we start 
from the equation 




dp 
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Gpm 
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(r 3 + > 



iU/3 
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(167) 



In this equation p is the pressure, m = m(r) is the 
matter mass inside a sphere of the radius r, p = p(r) is the 
matter density at the distance r from the center, r g is the 
function of m(r). 

Suppose the equation of state is p = Kp T , where K 
and T are constants. For numerical estimates we shall use 
their values from [p^]. 

For rough estimates we set p = Const and replace 
dp/dr by —p/R, where p is the average matter pressure and 
R is its radius. Under the circumstances we obtain from eq. 



p 

pc 2 



•g 
2R 



(168) 



If R ^> r g , then the term r g / f is negligible. Setting 
M « pR 3 , we find the mass of equilibrium states as a func- 
tion of p: 



M 



{K/Gfl^pP-WW. 



(169) 



It yields the maximum Chandrasekhar mass ^% M = 
(K/G) 3 ' 2 at p » p . 

However, according to eq. (168), there are also equilib- 
rium configurations at R < r g . In particular, at R -C r g 



we find from eq. (168) that the masses of the equilibrium 



configurations are given by 

M = c 9/2 1Q -l iir -3/4 Gf -3/2 /9 -(r-l/3)(3/4) > (17Q) 

These are the configurations with very large masses. For 
example, the following equilibrium configurations can be 
found: 



Figure 5: The plots of right-hand and left-hand sides of eq. 
(168) against mass M 



The following conclusions can be made after consider- 
ing the plots of the above kind: 

1 . For each value of p < p max there are two equilib- 
rium states (with R > r g and R < r g ). 

2. There are no equilibrium configurations whose den- 
sity is larger than a certain value p max ~ 10 16 g/cm 3 . (At 
the densities exceeding p max the curves do not intersect). 

More accurate conclusions about the internal structure 
of the configuration can be obtained from the equation of 
t he h ydrostatic equilibrium obtained for the gravity force 
(Ull): 
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where m = m(r) and r g 
is 

P 



i(r)). The equation of state 



rib 



P I c 2 , 



dn/dp ' 1 

where the baryons density rib as me function of p is given 
by the approximation Harrison equation pl[ ] which takes 
place from 7 to at least 10 16 g cm~ 3 : 



nb^Apil + Bp 1 ^)- 4 / 9 , 



(173) 



where A = 6.0228 ■ 10^ and B = 7.7483 ■ l(T iu in CGS 
units. 

In addition to the ordinary solution (i.e. configurations 
of the white dwarfs and neutron stars) there exist solutions 



i s 



with large masses. Fig. || shows an example of that kind 
of solutions for p(r). It is a configuration with the mass 
2.6 • 10 6 Af Q and the radius 0.378 • R . 



0.2 



0.4 0.6 0.8 

r/R 



Figure 6: The density distribution inside the object with the 
mass M = 3.9 ■ 10 35 .g and the radius R = 1.6 ■ 10 7 cm 

Are the configurations with large masses stable? The to- 
tal energy of the degenerate gas is E = Ei nt + E gr , where 
Ei n t is the intrinsic energy and E gr is the gravitational en- 
ergy. The gravitational energy of the sphere with the radius 



E gr = - I dm(r) x(r) m(r), 



(174) 



where 



X (r) = / dr> (r')- 2 (l - 1/7), 

J r 

r g = 2Gm(r)/c 2 , f = ( r g {rf + (r') 3 ) 1/3 , 



m(r) = 4tt / dr'p(r'Y{r'Y 
Jo 

The function x( r ) is given by 
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(175) 



(176) 



where T is the degenerated hypergeometric function. Ap- 
proximately 

X=(l/r)[(l-exp(-0.7r/r s )]. (177) 
Therefore, at p = Const up to a constant of the order one 

GM 2 



E — — - 



R 



[1 - exp(-0.7i?/r s )] 



(178) 



The intrinsic energy Ei nt — J u dm, where u is the 
energy per mass unit. For the used equation of state u = 
K(T — l) _1 p r_1 . Thus, up to the constant of the order of 
one 



r-i 



E = KM p 
GM 5/ V /3 [1 - exp(-QAf~ 2/ V 1/3 )], 



(179) 



where Q = 0.7c 2 /2G. As an example, fig^Jshows the plot 
of the function E = E(p) for the nonrelativistic degener- 
ated Fermi gas of the mass M = 2.5 • 10 38 g in comparison 
with the neutron star of the mass M = 10 33 g in fig. @. 
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Figure 7: The plot of the energy E vs. the density p for the 
configuration with the mass M = 2.5 • 10 6 M© 
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Figure 8: The plot of the energy E vs. the density p for the 
neutron star with mass 10 33 g 



The analysis of such plots shows that the function E = 
E(p) has the minimum. Thus, the above equilibrium states 
of the large masses are stable. 

More rigorous investigation that confirms this result 
was carried out in [?]. 



14 Conclusion 

It follows from the above results that the equations under 
consideration do not contradict available experimental data 
obtained in the Solar system. In paper | ^0| ] these equations 
were tested by the binary pulsar PSR19 13+16 and it was 
found out that the results are very close to the ones in gen- 
eral relativity. It is a consequence of the fact that the used 
distances from attracting masses are much larger than the 
Schwarzshild radius. At the conditions the function f(r) is 
very close to the radial distance r. However, the physical 
consequences between these equations are completely dif- 
ferent at the distances r <r g . The events horizon is absent. 
There can exist supermassive configurations of the degen- 
erated Fermi-gas. Candidates to the objects of such a kind 
are the galactic centers (IP4I -r- [38 1). 
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